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ABSTRACT 

The generation of metric coefficients of the coordinate 
transformation from a generally curved-sided domain boundary 
to the unit square (cube) is required for efficient solution 
algorithms in computational fluid mechanics. An algebraic 
procedure is presented for establishment of these data on the 
union of arbitrarily selected sub-domains of the global solution 
domain. A uniformly smooth progression of grid refinement is 
readily generated, including multiple specification of refined 
grids for a given macro-element domain discretization. The 
procedure is illustrated as generally applicable to non-simp1y 
connected domains in two- and three-dimensions. 

^Research principally supported by USAF Grant No. AFOSR-79-0005 . 
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COMPUTATIONAL REQUIREMENT 


Navier-Stokes Equations 


L(q,.) 

_ '^i 
at 

ax. 

Uj.qi + 

j 

= 0 

A ( q ^ ) 

= aiPi 

+ 32 

aq^ 

ax. j 

+ 33 = 0 


Coordinate Transformation 
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Numerical Solution Algorithm 
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Discussion 

The Navier-Stokes equations contain the vector divergence 
operator. The required transformation projects x. onto n.j with 
coordinate surfaces defined coincident with solution domain 
boundaries. The Cartesian description of dependent variables is 
retained, while the convection velocity is expressed in contra- 
variant scalar components. The numerical solution implementation 
requires nodal distributions of components of the forward and 
inverse Jacobians, and J, K, and are tensor summation indices. 


LOCALLY BI-QUADRATIC COORDINATE TRANS FORMAT I ON 


Physical Domain Transformed Domain 
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Three-Dimensional 

Discussion 
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The bi -quadratic cardinal basis {Nzin)} transforms the vertex 
and non-vertex node coordinate description of a smooth region of 
onto the unit square or cube spanned by the locally rectangular 
Cartesian coordinate system r[. The inverse transformation is 

non-singular for a range of non-midpoint definitions of the non- 
vertex node coordinates (x), yielding a non-uniform discretization 







DETAILS OF THE COORDINATE TRANSFORMATION 


Nodal Coordinates (XI}: 


Where; 
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Jacobians 
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Discussion 

Within a macro-domain, the components of both J and are 

continuous functions of nj and the global macro-node coordinate 
pairs (triples) {XI}, 1 £ I ^ n. Each global coordinate x^ pos- 
sesses an independent transformation; the corresponding Jacobian 
must be of rank n to assure existence of J”’^. Once the matrix 
elements of {XI} are defined, selection of any coordinate (m. na) 
defines a unique coordinate pair (xi , X 2 ), i.e., a mesh point on 
the refined grid in physical space. 
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